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SYNTHESIS  OF  TAYLOR  AND  BAYLISS 
PATTERNS  FOR  LINEAR  ANTENNA  ARRAYS 


INTRODUCTION 

The  requirement  for  low  sidelobes  from  array -type  antennas  is  a  long-standing  one.  The 
contributions  to  this  theory  extend  from  Dolph’s  utilization  of  Chebyshev  polynomials,  through 
Taylor’s  papers  on  linear  and  circular  apertures,  Bayliss’s  extension  to  difference-type  patterns,  and 
finally  to  recently  developed  techniques  which  provide  arbitrary  pattern  control  for  linear  arrays 
[1-8]. 


The  purpose  of  this  report  is  to  examine  some  of  the  more  recent  applications  of  these  syn¬ 
thesis  techniques  in  light  of  their  limitations  and  also  the  computational  capabilities  which  are  now 
available.  For  example,  at  the  time  Taylor  published  his  synthesis  procedure,  engineers  had  only 
slide  rules,  mathematical  tables,  and  mechanical  desk  calculators  to  generate  the  distribution  func¬ 
tions.  The  computational  capability  available  to  today’s  engineer  is  vastly  different,  and  we  will 
show  how  Taylor’s  and  Bayliss’s  procedures  can  be  modified  to  give  better  results. 

A  more  careful  look  at  the  synthesis  procedures  previously  mentoned  is  presented  in  Table  1. 

Dolph’s  synthesis  is  precise  and  gives  minimum  beamwidth  for  given  sidelobe  levels,  but  these 
constant  amplitude  sidelobes  are  not  desirable  for  larger  arrays  because  it  is  possible  to  radiate  most 
of  the  energy  into  the  sidelobes.  Taylor  solved  this  problem  by  allowing  the  far-out  sidelobes  to  fall 
off  as  dictated  by  an  amplitude  discontinuity  at  the  ends  of  the  aperture.  Taylor,  and  later  Bayliss, 
synthesized  continuous  distributions  and  sampled  these  to  obtain  array  excitations. 


Table  1  —  Synthesis  Procedures  for  Linear  Array  Apertures 


Procedure/ 

Date 

Continuous  or 
Discrete 

Limitations 

Dolph/47 

Discrete 

Poor  results  for  large  arrays 

Taylor/ 52 

Continuous 

Inexact  for  low  sidelobes,  small  arrays 

Bayliss/68 

Continuous 

Inexact  for  low  sidelobes,  small  arrays 

Hyneman/68 

Continuous 

Inexact  for  low  sidelobes,  small  arrays— iterative 

Stutzman/72 

Continuous 

Inexact  for  low  sidelobes,  small  arrays— iterative 

Elliott/76 

Continuous 

Inexact  for  low  sidelobes,  small  arrays— iterative 

Elliott/77 

Discrete 

Applies  all  continuous  procedures  to  discrete  arrays 

Manuscript  submitted  June  15,  1981. 
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Some  recent  applications  have  called  for  lower  sidelobes  and  smaller  arrays,  thereby  pressing 
the  limitations  of  the  Taylor  and  Bayliss  synthesis  procedures.  The  problem  of  discretizing  contin¬ 
uous  aperture  distributions  has  been  treated  [9-10] .  The  technique  used  in  this  report  is  different 
from  those  of  Winter  and  of  Elliott,  but  it  is  mathematically  related  to  Elliot’s  technique. 


REVIEW  OF  TAYLOR  SYNTHESIS  PROCEDURE 

A  brief  review  of  the  Taylor  synthesis  procedure  is  given  here.  The  key  to  this  procedure  is  the 
equal-sidelobe  pattern  function  which  is  the  continuous-aperture  analog  to  the  Chebyshev  polyno¬ 
mial  pattern  for  arrays: 


E(u)  =  COS  TTy/  u2  -  A2,  (1) 

where  u  =  7ra  sin  0 /X,  a  is  the  length  of  the  aperture  and  6  is  the  angle  measured  relative  to  the 
normal  to  the  array.  This  function  has  a  maximum  value  of  cosh  ttA  at  u  =  0  and  unit  sidelobes 
extending  to  u  =  ±  °°.  Taylor  showed  that  the  pattern  of  Eq.  (1)  is  not  physically  realizable  from  a 
continuous  aperture  distribution,  just  as  the  Dolph  array  excitation  becomes  increasingly  impracti¬ 
cal  in  the  limit  of  large  arrays.  His  brilliant  solution  to  this  problem  was: 

1.  For  all  zeros  of  the  synthesized  pattern  functions,  which  we  will  call  Es(u),  from  the  nth 
from  the  origin  to  °°,  the  locations  will  be  the  same  as  those  from  a  uniformly  illuminated  aperture 
of  the  same  size.  That  is, 

Eg{u)  =  0  for  u  =  n  for  n  >n. 

2.  For  the  first  7f  -  1  zeros,  their  locations  will  be  determined  by  the  zeros  of  E(u),  scaled  so 
that  the  nth  zero  is  located  at  u  -  7f. 

The  aperture  distribution  is  determined  by  performing  a  Woodward  synthesis  of  Es(u).  That  is, 
we  define  a  set  of  functions  of  the  form 

Fn(u)  =  sin  (u  -  n)ir/(u  -  n)n, 
and  then  construct  Eg(u )  from  the  Fn(u) 

oo 

Eg(u) =  F,(n)Fn(u).  (2) 

n  =  -°° 

Since  we  have  defined  £,(n)  =  0  for  n  >  n ,  Eq.  (3)  becomes 

7T  -  1 

£,<«)=  X)  EM)FnW-  (3) 

n“-S+l 
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Fourier  transformation  of  Eq.  (3)  yields  the  aperture  distribution: 


Es(u)ej2xunladu 


/  "  * 

5 


Es(n)Fn(u)ej2xu1,ladu. 


That  is,  A(x)  is  a  weighted  sum  of  integrals  of  the  form, 

(u  -  n)ir 


(4) 


Letting  u  =  u  -  n  results  in 


r  sin  u  tt 
J2nitx/a  I  — 

J „  UV 


ej2xu  nladu> 

Since  the  imaginary  part  of  the  integrand  is  odd,  this  becomes 

7r/a  du' 


ej2nitx /a  f  gin  cos  2xu’ 

I  u'lt 

OO 

nirx/a  f  1  [~sin  u'n(  1  -  2x/a)  +  sin  + 

J  2  L  uV 

OO 


2x/ a) 


du' . 


(5) 


A  standard  definite  integral  is 


=  7r  for  b  >  0 
=  0  for  b  =  0 
=  —i r  for  b  <  0 

Application  of  this  integral  to  Eq.  (5)  and  thence  to  Eq.  (4)  yields 

~n- 1 

A(x)  =  Eg(n)ej2nnxla 

n = - n+ 1 

n  - 1 

=  Eg(0)  +  2  £g(n)  cos  2irnx/a  for  |*|<a/2  (6) 

n  =  l 

=  0  for  | jc |  >  a/2. 


L 


sin  bzdz 


if n  “ 
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The  continuous  aperture  distribution  given  by  Eq.  (6)  is  sampled  to  give  the  element  excita¬ 
tion  values  for  a  discrete  array.  This  last  step  is  approximate,  and  the  pattern  function  of  the  array 
is  obviously  different  from  Es(u).  This  approximation  is  acceptable  provided  that  the  number  of 
elements  in  the  array  is  much  greater  than  n  and  the  sidelobe  level  is  not  extremely  low.  Figure  1 
is  an  example  of  a  case  in  which  the  synthesis  procedure  gives  an  unsatisfactory  result.  For  a  side- 
lobe  level  of  50  dB  below  mainbeam  and  n  =  8,  a  30-element  array  has  the  computed  pattern 
function  shown.  The  near-in  sidelobes  are  unduly  low,  whereas  the  first  eight  sidelobes  should 
be  about  the  same  level. 


ARRAY  PATTERN  FUNCTIONS  IN  TERMS  OF  ZEROS 

Elliott  used  a  synthesis  technique  which  relates  the  discrete  array  distribution  directly  with  the 
array  pattern  [9] .  We  also  use  this  relationship,  and  our  procedure  achieves  identical  results  with 
those  of  Elliott.  However,  the  actual  computations  are  different,  and  it  is  desirable  to  compare  the 
techniques. 

Elliott  expresses  the  pattern  function  as  a  polynomial  in  w,  where  w  =  2jr*/A)sin  6  The 
zeros  of  this  polynomial  are  given  by  wn,  which  are  normally  located  on  the  unit  circle.  Once  he 
has  the  wn  properly  adjusted,  he  completes  the  synthesis  by  multiplying  out  the  product  expression, 
n(u>  -  wn),  into  the  polynomial.  The  coefficients  of  the  polynomial  are  the  excitations  of  the 
array  elements. 


'-90.00  -60.00  -30.00  0.00  30.00  60.00 


0EG 

Fig.  1  —  Conventional  Taylor  synthesis,  N  =  30, 
3  «  8,  50-dB  sidelobes 
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Our  procedure  also  uses  the  pattern  function  zeros  in  a  product  expression.  Since  the 
patterns  are  symmetric,  our  expression  can  be  of  the  form,  II(cos  z  -  cos  zn ),  where 
z  =  (2tts/X)  sin  0.  We  cannot  multiply  this  product  expression  out  to  obtain  the  coefficients 
directly  since  we  require  terms  of  the  form  cos  nz  rather  than  cos  nz.  Rather,  we  carry  out  a 
synthesis  exactly  analogous  to  that  used  by  Taylor.  Uniformly  spaced  pattern  function  samples 
are  found  by  using  the  product  expression.  These  pattern  samples  are  used  in  a  Fourier  series  to 
find  the  array  illumination. 

The  procedure  relies  on  the  equivalent  location  of  pattern  function  zeros  for  the  line  source 
and  for  the  discrete  array.  Whereas  the  zeros  for  the  pattern  of  a  uniform  line  source  distribution  are 
located  at  u  =  n,  the  analogous  relationship  for  a  discrete  array  is  z  =  nir/N,  where  z  =  2jts  sin  0/A, 
where  s  is  element  spacing  and  N  is  the  number  of  elements  in  the  uniformly  excited  array. 

The  transformation  of  Taylor’s  procedure  is  easily  seen  to  consist  of  locating  the  zeros  in 
step  1  above  at  z  =  nir/N  for  n>H  and  then  scaling  the  first  n  zeros  of  Eq.  (1)  so  that  the  nth  zero 
is  located  at  z  =  nir/N. 

Appendix  A  lists  the  resulting  equations  for  Taylor  arrays  of  both  even  and  odd  N,  and 
Appendix  B  lists  the  equations  for  Bayliss  arrays  (yielding  monopulse  difference  patterns)  of  both 
even  and  odd  N.  Figure  2  is  an  example  of  a  Taylor  array  pattern  with  sidelobe  levels  of  50  dB  with 
7f  =  8  and  N  =  30.  These  equations  can  be  straightforwardly  programmed  for  automatic  processing 
by  a  digital  computer.  Many  programmable  calculators  now  have  sufficient  memory  to  implement 
these  programs.  Appendix  C  lists  programs  for  carrying  out  the  synthesis  and  evaluating  the 
pattern  functions  with  an  HP-41C  programmable  calculator. 


Fig.  2  —  Discretized  Taylor  synthesis, 
N  ■  30,  n  “  8,  50-dB  sidelobes 
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Appendix  A 

DESIGN  EQUATIONS  FOR  LINEAR  ARRAYS 
WITH  TAYLOR-TYPE  PATTERNS 


These  equations  will  determine  the  aperture  illumination  coefficients  for  a  linear  array  of 
N  elements  to  produce  a  Taylor-type  pattern  function  with  H  sidelobes  on  each  side  of  the  main 
beam  at  a  level  of  L  dB. 


This  design  procedure  involves  three  steps.  The  first  W  -  1  zeros  of  the  pattern  are  determined. 
Then  the  appropriate  pattern  function  samples  are  determined.  Finally,  the  array  element  illumina¬ 
tion  coefficients  are  determined  by  a  harmonic  analysis  of  the  pattern  function  samples. 

A  particular  advantage  of  this  synthesis  is  that  the  knowledge  of  all  of  the  pattern  function 
zeros  allows  the  computation  of  the  pattern  function  as  a  product  rather  than  as  a  polynomial.  The 
product  computation  involves  only  one  trigonometric  function  evaluation  for  each  pattern  func¬ 
tion  value.  All  other  constants  need  to  be  evaluated  only  once  for  each  array. 


The  pattern  function  zeros  are  given  by 


-  27rny/A2  +  (n-  1/2)2 
Ny/A2  +(n  -  1/2)2 


for  rt  =  1  to  n  -  1 


2irn 

N 


for  n  =  n  to  M  , 


(Ala) 

(Alb) 


where 

i1) 

and  A  is  given  by 

A  =  \  cosh- 1  [lO(L'20>J  (A2a) 

=  (L  +  6.02)/27.29,  (A2b) 

where  L  is  the  sidelobe  level  (positive)  in  dB.  Equation  ( A2b)  is  an  excellent  approximation,  espe¬ 
cially  for  large  L. 


M 


=  int  ^ 
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The  pattern  function  is  given  by 


E(z)  =  cos  - 
2 


M  /cos  2  -  COS  2  \ 

n(-T^r) 

M  /  COS  2  -  COS  2  \ 

n(-T^r) 


N  even 


Nodd 


(A3) 


The  pattern  samples  to  be  used  to  find  the  array  element  illumination  coefficients  are  given  by 

(2itm\ 


“m=£(- 


N  ) 


for  m  =  1  to  n  -  1  . 


( A4) 


The  element  excitation  coefficients  are  given  by 


ep  1+2  ”m 
m  =  1 


m(2p  -  l)ir 

cos -  Af  even,  p  =  1  to  Af  +  1 

N 


7T-  1 


+2Sa 


2  mptr 

cos - 

m  AT 


AT  odd,  p  =  0  to  M  , 


(A5) 


where  p  is  an  index  or  element  number  starting  at  the  center  and  moving  to  either  end  of  the  array. 
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DESIGN  EQUATIONS  FOR  LINEAR  ARRAYS  WITH 
BAYLISS-TYPE  DIFFERENCE  PATTERNS 


Appendix  A  gave  the  design  equations  for  linear  arrays  with  Taylor-type  patterns,  which 
produce  a  main  beam  with  slightly  larger  beamwidth  than  that  of  the  Dolph  synthesis  but  in  gen¬ 
eral  with  higher  gain.  In  some  applications;  such  as  monopulse,  we  might  require  a  difference  pat¬ 
tern.  Bayliss  presented  a  synthesis  procedure  for  difference  patterns,  analogous  to  that  of  Taylor. 
In  this  appendix  we  adapt  the  Bayliss  procedure  to  discrete  arrays. 

As  in  the  case  of  the  Taylor  synthesis,  the  application  of  discrete  arrays  involves  three  steps. 
The  first  7f  -  1  off-axis  zeros  of  the  pattern  are  determined.  Then  the  appropriate  pattern  function 
samples  are  determined.  Finally  the  array  element  illumination  coefficients  are  determined  by  a 
harmonic  analysis  of  the  pattern  function  samples. 

The  pattern  function  zeros  are  given  by 


27rM"  +  2 


= 


N~y A2  +  n  2 


2rr 


("+iV 


A2  +  n2 


Nyf A2  +  n2 

N 


for  n  -  1,  2,  3,  4 


for  n  =  5  to  n  -  1 


for  n  =  n  to  M 


(Bla) 


(Bib) 


(Blc) 


where 


M 


/  N  -  2  \ 

=  int(— ). 


In  this  case  it  is  necessary  to  find  both  A  and  qn  from  graphs  in  Bayliss’s  paper  [4] .  For  50  dB 
sidelobes,  A  =  2.42,  qj  =  2.78,  q2  =  3.18,  q3  =  3.85,  and  q4  =  4.65. 


The  pattern  function  is  given  by 


til  r  ,  21 A  r  *1  1 

E(z)  =  sin-  |  ]  I  cos  z- cos  z  I  /sin — Jj  cos — -cosz  I 
2  n  =  l  4  n  =  l  *-  2  -* 

M  Zj  M  T  Zj  I 

=  sin  z  rT  cos  z  -  cos  z  I  /sin — T1  cos — -cosz  I 
n  =  l  "j/  2  n  =  1  L  2  "J 


N  even 


N  odd 


(B2) 
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E{z)  is  normalized  to  unity  at  z  =  zx  / 3,  which  is  near  the  pattern  maximum.  If  a  more  precise 
pattern  maximum  is  desired,  a  better  multiplying  constant  can  easily  be  found. 

The  pattern  samples  to  be  used  to  find  the  array  element  illumination  coefficients  are 
given  by 


b 


m 


for  m  =  1  to  n  . 


The  element  excitation  coefficients  are  given  by 


(B3) 


„  ,  .  ff(2m  -  l)(2p  -  1) 

e  -  2  >  b  sin - 

P  ^  m  2 N 


for  N  even,  p  =  1  to  M  +  1 


=  2 


t 


b  sin 

m 


ir(2m  -  1  )p 
N 


for  N  odd,  p  =  1  to  M  +  1  (B4) 


where  p  is  an  index  of  the  element  number  starting  with  zero  at  the  center  of  the  array.  For  N  odd, 
the  center  element  of  the  array  always  has  zero  excitation.  The  excitations  on  one  side  of  the  array 
are  the  negative  of  those  on  the  other  side. 
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PROGRAMS  FOR  THE  HP-41C  CALCULATOR 


This  appendix  presents  programs  for  the  HP-41C  calculator  for  the  design  equations  of  Appen¬ 
dices  A  and  B.  The  software  consists  of  four  programs,  SUM,  DIF,  IN.  and  SL.  “SUM”  contains  the 
equations  for  synthesizing  Taylor-type  sum  patterns;  “DIF"  contains  equations  for  Bayliss-typc 
difference  patterns;  “IN”  contains  subroutines  that  are  used  by  both  programs;  and  “SL”  is  a 
routine  for  calculating  the  peaks  of  the  sidelobes  of  the  synthesized  array.  The  number  of  registers 
used  by  the  programs  and  the  number  of  card  sides  required  for  storage  are: 


Program 

Registers 

Card  Sides 

SUM 

30 

2 

DIF 

42 

3 

IN 

39 

3 

SL 

19 

2 

130 

10  (5  cards 

It  is  possible  to  synthesize  aperture  distributions  using  either  SUM  and  IN  or  DIF  and  IN.  These 
programs  require  at  least  one  additional  memory  module.  Furthermore,  the  programs  use  nine 
registers  for  variables,  indices,  and  constants.  Table  Cl  correlates  the  number  of  registers  avail¬ 
able  for  synthesis  parameters  with  the  number  of  additional  memory  modules  in  use.  The  avail¬ 
able  registers  are  used  for  the  pattern  samples  am  and  bm  and  for  the  pattern  function  zeros 
(cosines)  zp.  The  number  of  these  registers  is  7T  +  M.  Therefore,  the  size  of  array  that  can  be 
synthesized  for  any  given  configuration  of  Table  Cl  depends  on  n  .  For  a  50-dB  sidelobe  require¬ 
ment,  7f  will  be  about  8.  Roughly  speaking,  an  array  of  55  to  65  elements  for  difference  and  80  to 
90  for  sum  can  be  synthesized  using  one  memory  module  by  trading  programs  in  and  out  of  the 
machine,  and  an  array  of  90  to  100  elements  can  be  synthesized  with  all  programs  loaded  using 
two  modules.  The  maximum  array  size  that  can  be  handled  using  three  modules  is  310  to  320  for 
difference  and  about  340  for  sum.  It  appears  that  one  or  two  memory  modules  should  suffice  for 
most  requirements. 


Table  Cl  —  Registers  Available  after  Loading 
Indicated  Program  Complements 


Program  Complement 

Number  of  Memory  Modules 

1 

2 

3 

SUM  +  IN 

48 

112 

176 

DIF  +  IN 

35 

99 

163 

SUM  +  DIF  +  IN 

11 

71 

135 

SUM  +  DIF  +  IN  +  SL 

— 

54 

118 
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The  procedure  for  running  the  programs  is: 

1.  Allocate  memory  by  XEQ  SIZE  (9  +  M  +  n). 

2.  Load  the  appropriate  program  complement. 

3.  Enter  either  XEQ  SUM  or  XEQ  DIF. 

4.  The  display  will  prompt  for  N,  L,  and  NBAR.  N  can  be  even  or  odd.  L  is  sidelobe  level  in 
positive  dB.  DIF  will  also  prompt  for  A,  Ql,  Q2,  Q3,  Q4. 

5.  After  calculating  zn  and  loading  cos  zn  into  registers  starting  with  (9  +  n),  the  display  will  ask 
whether  you  want  a  listing  of  peak  sidelobes  (SL)  or  aperture  distribution  (EP).  After  the  sidelobes 
or  excitation  coefficients  are  listed,  the  display  will  ask  whether  you  want  the  other  set  of  parameters 
calculated  and  listed. 

The  routine  SL  computes  the  sidelobe  level  relative  to  the  main  beam  level  by  evaluating  the 
pattern  value  at  a  point  midway  between  pattern  zeros.  This  computation  is  admittedly  approximate 
because  the  pattern  maximum  is  in  general  not  exactly  midway  between  zeros.  The  main  beam  pat¬ 
tern  value  is  computed  for  z  =  0.  The  difference  pattern  maximum  is  computed  for  z  =  Zj  / 3.  This 
factor  was  found  to  be  accurate  for  50  dB  sidelobes.  The  exact  multiplying  factor  will  be  somewhat 
larger  for  higher  sidelobes  (L  <  50),  and  it  can  be  found  quickly  by  obtaining  z1  and  executing  PA: 

RCL  (9  +  n)  gives  cos  z1 
ACOS  gives  z1 

k  new  multiplying  factor,  such  as  .4 

* 

COS 
STO02 
XEQ  PA  . 

Alternatively,  k  can  be  found  from  Fig.  4  of  Baylissci ,  which  defines  the  beam  maximum  by  pQ, 
where  k  -  pQ  /  §  1 .  ( §  x  corresponds  to  our  z  1 ) 

Once  the  desired  value  of  k  has  been  found,  go  to  lines  110,  111  in  DIF,  and  exchange  ft, 

*  for  3,/.  It  is  now  necessary  to  reload  the  reference  main  beam  pattern  value  into  R08.  This  calcu¬ 
lation  starts  at  line  61  of  SUM  and  105  of  DIF.  Alternatively,  you  can  simply  rerun  the  program. 

The  pattern  value,  in  voltage  and  normalized  to  mainbeam  level,  is  found  by  keying  in  the 
value  of  z  in  degrees,  then  keying  COS,  STO  02,  XEQ  PA. 

The  registers  used  are: 

00  N 

01  M 

02  A2  and  cos  zm  for  PA 

03  n 

04,  05  loop  indices 

06  multiplying  constants 
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07  accumulator  for  E(z),  ep 

08  main  beam  reference  value 

09  to  (8  +  n)  computed  values  of  am ,  bm 

(9  +  n)  to  {8  +  n  +  M)  computed  values  of  cos  zn. 

Program  IN  contains  the  following  subroutines: 

IN  Asks  for  input  data  N,  L,  NBAR 

ZN  Completes  calculation  and  storage  of  cos  zn 

BR  Asks  for  choice  of  sidelobes  or  aperture  distribution  and  branches 

to  EP  or  SL 

PR  Prints  element  excitations  ep 

PA  Computes  pattern  value  for  am,  bm,  or  SL  routines 

EP  Completes  calculation  of  ep. 

The  programs  use  flags  00  and  01  to  indicate  the  following  conditions: 

Flag  00  is  set  for  N  even 
clear  for  N  odd 

Flag  01  is  set  for  DIF  execution 

clear  for  SUM  execution. 

The  use  of  registers  by  program  PA  precludes  the  use  of  the  plot  subroutines  resident  in  the 
printer. 

Note  that  the  sidelobes  and  pattern  values  obtained  with  these  programs  are  all  relative  to  the 
main  beam  level.  No  information  concerning  gain  or  aperture  illumination  efficiency  is  computed. 
The  aperture  distribution  can  be  used  to  compute  aperture  efficiency  or  gain . 

The  programs  and  sample  printouts  are  listed  on  the  following  pages. 
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